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1. Introduction 


The double charge exchange (DCX) reaction with pions is still attracting a good deal 
of attention (see ref. [1] the latest review of the subject). There was hope that special 


information about an affection of one nucleus by another would appear in the reaction. 


- Unfortunately, the meson factories were built in the 1970's, any type of such correlations 


was masked by other complex interactions at highly accessible meson energies at that 


time. Only recently [2-4] as small as 50 MeV meson energies are available and nucleon 


correlations have begun to surface. 
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The double charge exchange process necessarily involves two neutrons. The ingoing 
pion r* has successive charge-exchange on a neutron. In a second step the produced pion 
x? exchanges the charge to another neutron to produce a proton and a pion x~, which then 
leaves the nucleus and is detected. Theorists expect in general that the pions interact in the 
above picture only if the DCX process takes place on valence nucleons which are within 1 fm 
(about nucleon radius ) of each other. In other words the strong short-range correlations to 
nuclear wave functions mist be involved if one wants to reduce a considerable disagreement 
with the experiment [2]. In fact the theories including correlation effects come much closer 
to the data and a fairly satisfactory agreement is reported [3-4]. Another regularity of the 
reaction is the sensitivity to the role played by antisymmetry of nuclear wave functions, 
which allows the nucleons of the core to participate in an active way in the reaction [5]. - 
The alleged role played by quarks in the DCX reaction has added elements of interest to | 
the field [6]. 

In all reported calculations the strong limited shell model or even uncorrelated free 
nucleon wave functions were used. This limitation imposes applications to the light nuclei 
and/or to nuclei with a few nucleons out of the closed shell only. Bleszynski and Glauber’s 
wave function of the '4C nucleus is assumed to separate into a wave function for two 
valence neutrons in the p-shell orbit and a wave function for the core. The authors argue 
that the core plays a passive role and they omit its wave function in derivations. The 
two p-shell nucleons, of the total angular momentum zero, can only be in one of the two 
states, the spin singlet 1S or the spin triplet *P, state and final nuclear wave function 
is a free combination of these two states. For description of the calcium isotopes the 
same authors have used pure f7/ą wave functions with the parametrized p-configuration 
admixture. More fundamental approach, the seniority scheme was proposed for description 
of the DCX reaction by Auerbach, Gibbs and Ginocchio [4], but in principle the model is 
only valid for nuclei with a few valence particles or holes of one kind, for example, for the 
calcium isotopes. | 

In the present paper we apply the random phase approximation based on the togi 
particles formalism (QRPA) to give a consistent approach to the DCX reaction in a wide 


mass region of nuclei. The motivation for doing results from the observation that the 
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QRPA formalism provides a very useful way of description of the beta and double beta 
decay processes [7-9]. Because in the DCX reaction the nuclear structure is as important 
as in the beta decay, one can expect similar usefulness of the QRPA. Moreover, further 
formal analogies between operators cause that both reactions show some similarities and 
have already been discussed in literature [10]. 

The article is arranged as follows. A form of the DCX operator for the p-wave pion in 
non-static and static limit is studied in Sec. 2. The basic features of the QRPA approach are 
reviewed in Sec. 3. Next section gives formulae for the DCX total amplitude in the ground 
state to ground state transition case. Summary and remarks (Sec. 5) and appendices close 


the paper. 
2. The Double Charge Exchange Operator 


We will consider the simplest local xNN interaction Lagrangian linear in the pseudo- 
scalar-isovector pion field p. The Lagrangian is a scalar-isoscalar quantity, the nucleon 
fields must therefore appear in a pseudoscalar-isovector combination. This leads to the 
well-known pseudoscalar coupling Lagrangian 


‘Lps(z) = -gd(zjiqy,FY(z)-6 . (1) 


The experimentał value of the coupling constant g has been determined accurateły from 
pion-nucleon and nucleon-nucleon scattering g? = 180. The Lorentz-invariant form of the 


„xNN interaction deduced from the Lagrangian (1) is (e.g. [11]) 


(N (P')|CPsiN (p) ,x* (k)) =igvŻi(p)Fyu(p): 6 . (2) 
In the above equations a tilde indicates an isovector character of the operators. Since 
we deal with nucleons which are not far off-shell, we use the standard expression for the 
nucleon spinors u to obtain 


"e _ (E+M E+M)? | op BR. 
eeu) = (FAM ele - ale © 
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Here € is the Pauli spinor for the nucleon, M - the nucleon mass and energy transfer 
w = E' — E. The total energy of nucleon fulfills the usual relativistic relation E = E, = 
(p? + M?)2. An analogous expression one can write for the quantities E' and p'. If we 


expand eq. (3) to the lowest order in (p’ — p) and (p' + p) we find the result 


W (P) eas), +4 N) = ize o- ASF EE | (4a) 


= Le. (F-n) EA (48) 
- The second line involves a new interaction constant originated from the second form of the 
«NN interaction, the pseudovector coupling (PV)!. Its experimentally determined value 
is f? m 1.005. The operator Vn = 1 (Fx ~ Fn) acts to the right and left on a single 
nucleon and p' — p = k is transfer of pion momentum. 

The static non-relativistic reduction of the Lps Lagrangian is easily shown to lead to 


the effecti ve pion-nucleon Hamiltonian 
H f (o ) G P ) p (5) 
ei N m z 


'The Hamiltonian (5) is linear in the pion momentum, and therefore it generatęs p-wave 
pion nucleon interaction only?. 

Egs. (4) and (5) need some comments. At the most fundamental level the whole picture 
of the DCX reaction depends on the current choice of the NN interaction operator. In 
the earliest calculations and many of the recent ones the simple static operator (5) has 
been used. This operator was obtained by keeping just the leading term in the non- 
relativistic limit of the pseudoscalar coupling (1), which is conventionally used to describe. 
the r NN vertex. If we deal with some relativistic corrections introduced by some particular 
approximations, they can involve different correction terms into the original expression (3). 
There are at least three ways to switch on A pproziińzikida : (i), the Foldy-Wouthuysen 


transformation, (ii), Galilean invariance of the non-relativistic expression and (iii), the 


1 One can show equivalence of the both types of couplings in the non-relativistic limit and we will not 
pay attention to the PV form of the interaction. For details see e.g. [11]. 

2 Possible s-wave r — N interaction is based on mechanisms completely different from those described 
above and we do not intend to discuss this point in the present paper. 
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higher-order non-relativistic reduction of eq. (3). The final forms differ by some factors 
{13]. Furthermore, the whole picture must be AA as the basic process r N N requires 
an off-shell leg and so must be embedded in some more complicated interactions. All 
these ambiguities suggest a need for a more consistent relativistic picture, but because 
of difficulties in describing nuclei relativistically, we limit our considerations to the non- 
relativistic static Hamiltonian (5) only. 

A description of rNN scattering must correctly incorporate the contributions of se- 
cond-order processes generated by absorption or emission of a pion on the nucleon. One 
has to calculate amplitudes of these processes and therefore we need the matrix elements 


connecting the vacuum and a one-pion state |ra(q)) (e.g. [11]) : 


(r0(@!) Eunn) = iL ogne : (6a) 
(Olenwira(g)) = ~iLo-aree'™ x (6b) 


The transition matrix T for the direct and crossed terms shown in Fig. 1 is obtained using 


the standard second-order perturbation theory with the operator HzNN. 


3. Quasiparticle Random Phase Approximation and the form of the Hann 


Hamiltonian i 


The proton-neutron quasiparticle RPA formalism has been extensively discussed in 
literature [14-16], so we need only to present its basic feature here. 

The pn-QRPA was developed for the description of charge-changing excitations, such as 
the Gamow-Teller transitions or charge exchange with pions. In these processes iransitions 
from the 0* ground state of an even-even nucleus (Z,N) to states in the odd-odd nucleus 
(Z+1,N#1) are considered. The QRPA phonons with angular momentum J and its third 
component M are defined by the ansatz (e.g. [17]): 


Qh = E [XG (PnJM)-Yn,Ct(pnJM)] , (7) 


pn 
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Figure 1. Tlustration of the direct (a) and crossed (b) contributions. 


where indices p and n distinguish between proton and neutron states. The proton-neutron 
pair-creation and annihilation operators are defined with the aid of the Clebsch-Gordan 


coefficients 


Ct(pnJM) =)" (jgmpinmnlIM abn bhm, (8) 


Mp: My 


and C(poJM)=[Ct(pnJM)]'. The forward- and backward-going amplitudes Xia) and 
Ym) z are obtained by solving the QPRA eigenproblem 


MIE Pie = o 


2 (cs 1) 2 (ra) sę (10) 


The submatrices A and B are expressed explicitly 


with normalization 


AJM iq! = (BCSIC(pnJM)H,„C'(p'n'J'M')|BCS) = 
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= bpp! inn OOM M! (Ep tEn) iF 
—2gpp((Pn)J M|G\(p'n')JM)(upunup un + VpUndp Vn’ )— 


—2gn((pń”*)JM|GI(p'n'"")J My(uzponupow + ypintpuw) , (11) 


B/M = (BCSIC(pnJM)Ć(p'n' JM)H,,|BCS) = 


= —29pp((pn)J M |GI(p'n')J M)(upuntp Un + Vpyntuprun')- 
—2gpn((pn"")JMIGI(p'n" ")JM)(upyntprun: +UpUnUpYw) . (12) 


Here H,, is the nuclear Hamiltonian after the Valatin-Bogoliubov quasiparticle transfor- 
mation. The quasiparticle energies e, (a= p or n) and the occupation amplitudes ug, Va . 
are obtained by solving the BCS equations in some modeł space of a single-particle po- 
tential, the Wood-Saxon potential for example. The matrix elements of particle-particle 
((ab) J M|G|(a’t!)JM) and particle-hole ((ab71)JM|G|(a'b'"') interaction are related to 
each other through the Pandya transformation and can be obtained from realistic or phe- 


nomenological forces? 


. The constants gpp and gpa parametrize the two-body matrix el- 
ements ( |G| ) in order to take into account the finite model space restriction always 
involved into calculations and a renormalization for finiteness of a nucleus. The submatrix 
A in eq. (11) is responsible for the mixing of various particle-particle states in the odd-odd 
nucleus reached by the charge-changing transition because its elements are dominated by 
the particle-hole interaction and quasiparticle energies. The contribution of the particle- 
particle interaction is smaller implying that this interaction has only a little effect on the 
mixing of particle-hole states in the (Z+1,N-1) intermediate nucleus. On the other hand, 
the elements of the submatrix B (eq. (12) consist of the both interaction terms of the 
same order. This indicates that the particle-particle interaction can significantly enhance 
the ground state correlations by increasing the backward-going amplitudes Ye) z {19}. 


3 A realistic one-boson exchange potential was applied in papers {8, 9, 18] and a zero-range spin-isospin 
force in ref. [7], for example. 
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The supermatrix in eq. (9) has the dimension 2r, where r is the number od different 
possible proton-proton two-quasiparticle states of the angular momentum J and parity 
x(a = +). The diagonalization of this matrix gives r eigenvectors for each J* combination 
corresponding to the r positive eigenvalues Q, (1=1,2,...,r) and the same number negative 
eigenvalues, which we treat as nonphysical ones. The physical eigenvectors 


|mJ*M) = QT}, IRPA) (13) 


with the RPA‘ phonon operator QT given by eq. (7), can be constructed for two separate 
sets of states of the same intermediate nucleus (Z+1,N-1) depending on if the QRPA 
calculation is based on the parent (Z,N) or (Z+2,N-2) daughter nucleus ground states (see 
Fig. 2). 

Both sets of intermediate states should be physically identical. In fact they result from 
the two different QRPA calculations which causes their mathematical nonequivalence. 
Particularly they fulfill the orthogonality relation approximately in the sense 


(mJm|m' Im) 30 for |EX —EJ,| > AB 


and 


Fa [Teal Fons) [? ml 
JEZ EJ „| <AE 


with AE being of the order 1 MeV [19]. 

'To describe the DCX reaction in the pn-QRPA model one has to expand the tran- 
sition operator (5) in the quasiparticle operators ala, (bi ba) . We get (for details see 
Appendix A) 

p af + 
ky(k) = —Vi2— D ((-17Yju()|DOGy(KRZA || (14): 
JM pn 


Here the QRPA part of the transition operator nz can be expressed in the form: 
RIM = uptnCl(pnJM) + vpttnC(pnJM) + upunD(pnJM) — vpvnDt(pnJM) . (15) 


The Ct(pnJM) and C(pnJM) operators have been defined in eqs. (8a)-(8b) and other 


two operators are expressed as follows 


D(pnJM) = YigmpjnmalTM aby (—1)**™ aang, 
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(Z,N) (Z+1,N-1) (Z+2,N-2) 


(a) i (b) 


Figure 2. The QRPA model for the DCX reaction. One starts with the RPA proton-particle neutron- 
hole (a) or neutron-particle proton-hole (b) calculation to the intermediate (Z+1, N-1) nucleus states. 
The mode (a) corresponds to the ground state of the (Z,N) nucleus while mode (b) to the ground state of 
the (Z+2, N-2) nucleus. A transition from an unoccupied to an occupied single particle state is strongly 
stifled. 
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. D(pnJM) = [Dł(prJM)]' . (16) 


The tilde operation A(abJ M) = (—1)7+M A(abJM) is used here to have for the annihi- 
lation operator the usual phase convention of spherical tensors. Yym((Qk) is the spherical 
harmonic of degree J {25] and an expression for the GZM (k) form-factor can be found in 
Appendix A (eq. A6). 


4. The QRPA Model for the DCX Reaction 


Treating the DCX reaction we assume a sequential mechanism. In this picture the 
incident pion w+ of the beam exchanges its charge to a neutron. In the second step 
the produced pion x” exchanges the charge to another neutron to produce a proton and 
a pion m”, which then leaves the nucleus and is detected. During this process spin- 
isospin dependent nuclear interactions mix the unperturbed particle-hole states giving a 
contribution to meeieaa ese pair excitations (Fig.3). In the RPA scheme the basic DCX 
process can be described by the two-fold action of the one-exchange QRPA operator h, 
(14). Schematic picture of the mechanism is shown in Fig. 4. 

ln the second-order perturbation theory we can express the DCX transition amplitude 


Fike) JU Lx Vil Ml! 


mm’ 


+ crossed term. i (17) 


jt,7*(k)) and |f,r"(K')) are the ground states of the parent and daughter nuclei; k 
and k’ are in- and outgoing pion momenta and ez = (k? + +m2)? is the incident pion 
energy. Note that the summation in eq. (17) involves products of the transition matrix 
elements of the both QRPA mode calculations (compare Fig. 2 and remarks in Sec. 3). 
This trick allows to write eq. (17) symmetrically in both indices m and m’ because of 
quasi-orthogonality of the sets {|mJM)} and f|m'JM)). The procedure introduces an 


uncertainty in the energy denominator for which we involved an average value of the 


DCX Ground State Transitions... 141 


x 
z 
otites” oTe'** 
> zj Z 
Z a 
hole N + w hole N ma 
/ Bte't* ska "7 
- Z 
x 
(a) (b) 


Figure 3. Diagrammatic representation of the nucleon-hole contribution in the lowest order: (a) — direct 
term, (b) — crossed term. 


s fe 
QRPA excitation energies EZ and EJ, however this uncertainty is small compared to the 
mean value of the denominator. 

Using eq. (6) and the explicit form of the hp operator (eqs. (14)-(16) and (A3)) after 


cumbersome derivations one can find the final formula for the total amplitude F'* 


F(k,k') = $ Fy(k,k') ; (18a) 
J 
; it z C= J|m'J 
Pi(k,k’) = -> (4 ) Palcos ba) > a OC. 
mm’ ie. ng a Ró 


* See details in Appendix B 
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n”! p 


Figure 4. The DCX mechanism in the RPA model. Intermediate ststes of the (Z+1, N-1) nucleus are 
described in Sec. 3. > 
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; { [v > GEJ" (Ega Pylia + Yn) 8,51 ) | 
pn 


. [vi F Gia 2) (Xs. sipte + YGosen)] | . (18b) 
pn ‘ 


Here Py(cos 6,4) is the Legendre polynomial of the J* order [20] and the overlap between 
two J” states belonging to two different sets, (mJM|m’JM) is given by 


(mIM|m'JM) = Y (Romy Xian; — Fzwystiesya) OF OF) - (19) 


pn 


In eq. (19) the factor (0+|0+) when approximated by the BCS states is 


(0F oy") = CZ P tpp) | [(untn +Undn) - (20) 
p n 
In the above expressions all quantities with a bar correspond to the daughter (Z+2,N-2) 
nucleus and without a bar to the parent (Z,N) nucleus. 

The crossed term (Fig. 1b) is treated similarly. The corresponding amplitude satisfies 
crossing invariance: the crossed term is obtained from the direct term by interchanging 
the in- and outgoing pion momenta as well as energies: q +— —q'; Ex +— —tw. 

The scattering amplitude F(k, k’) is connected directly with the differential cross sec- 


tion 


(w) =|FlGKJE | (21) 


The accuracy of the single and double charge exchange data requires the careful treat- 
ment of kinematic factors. The differential cross section (21) must be evaluated in the pion- 
nucleus centre-of-mass frame, which approximately coincides with the laboratory frame. 
The amplitudes (18) are defined in the tN centre-of-mass frame. It is therefore necessary 
to incorporate the kinematic factors relating these two frames of reference. This can be 
done in a manner of any textbook (e.g. [11]) and results in a slight change of angular 
momenta k, k’ and an angle Są: [11,23]. 
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5. Summary and Final Remarks 


In the paper we have introduced formulae for the DCX transition amplitude within the 
pn-QRPA formalism. The static limit for the p-wave pion-nucleon effective interaction was 
used but in principle there are no difficulties to extend this approach onto relativistically 
corrected expression (4). Existing suggestions [23] and the author's preliminary estimations 
made for calcium isotopes [24] allow to admit that the recoil term is of no quantitative 
importance for p-wave pion-nucleon interaction. 

An extension of the DCX description which one has to involve in the low pion energy 
region is an inclusion of the s-wave pion effective interaction [11, 22]. Suitable formulae 
for the double isobaric analog transitions can be found elsewhere [24] and for the ground 
state transitions in a subsequent paper. 

The author would like to thank Professor Amand Faessler for suggesting the interesting 
problem of the DCX reaction with pions and for many helpful discussions. The fruitful 
advice and stimulating discussions on the QRPA formalism with Professor $.B.Khadkikar 
and Dr. J.Suhonen are also OPIN 

The work was supported by the Bundesministerium fiir Forschung und Technologie. 


Appendix A. 


We need to express the transition operator (5) in the second quantization formalism. 


Firstly using the relation (6) one can evaluate the pion operator part to have 
(ra) = (0lEnn(z)fr"(q)) = —VBi Fog cry, (41) 


where ry. is the nuclear isospin raising operator and we used pion plain wave expansion 
for the pion field operator ø. Next we define the transition operator hy from its density 
hp = f d?zW'(x)(p(x)W(x) and express it in terms of the creation and annihilation 


nucleon operators cf,cp (ch,cn) as 


h= -Vii LE || decPovsix,Qo-ae*Hx,0)| chen 5 > (42) 
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or after transformation to the quasiparticles 


h(a =-V2i LD FAM (gem |. (43) 
* on, JM 
The quasiparticle operator RJM is this of eq.(15) and FJM (q) is given by 
FAM (a) = DV "etiymyjn=malTM | | Bzdgyzta. Oradea]. (44) 
(m) | 


In the above equations € means additional intrinsic spin coordinates. Cumbersome calcu- 


lations allow to find more useful expression for F3 (q): 


FIM (q) = Vår V6 Yju(%q)GFn(9) > (45) 
where : 
J : zara a te tt, =t" , $ lp Jp | 
Gja = (-1)**3,L jain 3 (-1) s (1,04,0|:"0) R), (4) 7 r K | . (46) 
pr 1 r | 


Here Ys; is the spherical harmonic [25] and ja = /2ja +1. The symbols ( | ) and 
£ } mean the Clebsch-Gordan coefficient and 6j symbol, respectively. The formulae for 
the RE (a) function can differ for different radial dependence of nucleon wave functions 
p(x, €)(n(x, £)). In the case of harmonic oscillator wave functions one get the analytical 


expression 
healt " Hol, pła=n,=n ( qe 
= Ynp!inn!(2n, + 2, + 1)!(2nn + ln + L)!!2irt"=+=">="= g exp ( CH 

Ry Re (—2)er# (2, i li = t") +s+ s')! 


s ee Nt U kdo?) N, "+4 
s!(np — s)!'(2ł, + 2s + 1)!!s'!(n, — s')!(2/, + 2s! + 1)! 


Ls * (bq”j4) | (AT) 


s,8'=0 


The constant b is the harmonic oscillator length b = (h/Mw)3 = 1.0064 fm and Lg 
means the Laguerre polynomial [20] with 8 = }(lp + ln- Ù) +s +5’. 
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To express the amplitude F(k, k’) (eq.(6)) in terms of the QRPA quantities one needs 
to calculate two matrix elements (f,x~|h,|mJM) and (m'JM|hp|i,x*). It can be done 


by using the quasiboson approximation [26] 
(RPA| [C(abJM),C'(a'b'J'M')] |RPA) = 6.ar0wr65. p óm M! (B1) 


and the fact that the QRPA phonon annihilation operator Q7 gives zero operating on 


the RPA vacuum. Finally one can receive 
(RPA| [REM CtpnJ M) [RPA) = (-1)7+M ypunbpprónnójnóm=m „+  (B2) 


(RPA| [RM C(pnsM)| |RPA) =-(-1)tMupnbzpónnójróm=m -  (B3) 


With the aid of the formulae (B1)-(B3) one finds expressions for both needed matrix 


elements 


(f,9*;z7(k')|h„|mJM) = -v12 LV mO) 


` [Ee G;„(k') (Xo) s2ptin + Yer.) repo ) (B4) 
pn u 
and 
(en! JM|hgli, OF; (k) = VB 1 VE- Yanu) 
 |ZGzlk) (X spre + YM) spun) (B5): 
. pn «l 


The quantities with and without bars correspond to the daughter and parent nucleus 


ground state, respectively. 


Abstract: The quasiparticle random phase approximation (QRPA) was applied to.the 
double charge exchange (DCX) reaction with pions r* + zXN — #7 + Z+2XN-2. A 


case of the ground state to ground state transition is discussed. Final formulae for the 
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total amplitude and differential cross section are given within the p-wave r — N effective 


interaction. 


Streszczenie: W artykule zastosowano formalizm RPA do opisu reakcji z pionami podwój- 

nej wymiamy ładunku na jądrach atomowych. Dyskutowany jest kanał reakcji, w którym 

końcowe jądro atomowe znajduje się w stanie podstawowym. Otrzymano wyrażenia na 

całkowitą amplitudę reakcji i różniczkowe przekroje czynne dla efektywnego oddziaływania 

pion-nukleon w przybliżeniu fali p. 
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